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Abstract 

The Filament Based Lamellipodium Model (FBLM) is a two-phase 
two-dimensional continuum model, describing the dynamcis of two in¬ 
teracting families of locally parallel actin filaments EH- It contains 
accounts of the filaments’ bending stiffness, of adhesion to the sub¬ 
strate, and of cross-links connecting the two families. 

An extension of the model is presented with contributions from 
nucleation of filaments by branching, from capping, from contraction 
by actin-myosin interaction, and from a pressure-like repulsion between 
parallel filaments due to Coulomb interaction. The effect of a chemoat¬ 
tractant is described by a simple signal transduction model influencing 
the polymerization speed. Simulations with the extended model show 
its potential for describing various moving cell shapes, depending on 
the signal transduction procedure, and for predicting transients be¬ 
tween nonmoving and moving states as well as changes of direction. 
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1 Introduction 


In contact to flat adhesive substrates, many cell types tend to develop thin 
protrusions, called lamellipodia [35]. These are very dynamic structures, 
supported by a network of filaments of polymerized actin, which is continu¬ 
ously remodelled by polymerization and depolymerization as well as by the 
building and breaking of cross-links and adhesive bonds to the substrate. Po¬ 
larization of cells due to internal instabilities im inn] or to external signals 
[57 ] |8[ [12, 57] might lead to crawling movement along the substrate. This 
movement serves as a model for the motility of various cell types in natural 
environments, such as fibroblasts, tumor cells, leukocytes, keratocytes, and 
others. 

The dynamics of the filament network is a complicated process, and ef¬ 
fects to be taken into account, additionally to what has been mentioned 
above, include the nucleation of new filaments by branching off existing fila¬ 
ments, deactivation of filaments by capping, and contraction by actin-myosin 
interaction ([16]. references therein, and more references below). Various 
attempts have dealt with the formulation of mechanical and, consequen¬ 
tially, mathematical models for the involved subprocesses as well as for the 
whole integrated system ( [24] and references therein). On the level of in¬ 
dividual actin filaments, polymerization, depolymerization, branching, and 
capping are typically modeled as stochastic processes, where the regulation 
of polymerization as the key process pushing the lamellipodium outward 
has received the biggest attention |25[ 155] . Models based on individual fil¬ 
aments have provided possible explanations for various phenomena, such 
as the motility of pathogens in host cells m- For the description of the 
morphology dynamics of whole lamellipodia these models are too complex, 
however. Therefore, continuum models for the mechanics of the filament net¬ 
work have been used, where the choice of model is typically guided by the 
expected rheological properties, such as viscoelasticity or active contraction 
due to actin-myosin interaction mmmm- 

This work is a continuation of previous efforts |31j to systematically de¬ 
rive a continuum model from filament based descriptions by an averaging 
process similar to homogenization of materials with microstructure. This 
allows to include detailed knowledge or assumptions on all subprocesses. 
We discuss the modeling assumptions of the Filament Based Lamellipodium 
Model (FBLM), starting with those aspects, which are taken from |31| with¬ 
out changes. 

Geometry: As common in homogenization, the averaging process is fa¬ 
cilitated by idealizing assumptions on the microstructure. For lamellipodia, 
the two main assumptions are a restriction to a two-dimensional model mo¬ 
tivated by the observed small aspect ratios (100-200 nm thickness, tens of 
microns lateral extension, [35]) and the idealization to a network consisting 
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of two families of locally parallel filaments crossing each other transversally. 
The latter assumption is supported by experimental evidence for steadily 
moving cells |55| . It has to be conceded, however, that it is questionable in 
certain conditions such as retracting lamellipodia |15| . 

It is assumed that the whole cell is surrounded by the lamellipodium, 
whose width might vary along the cell periphery. Mathematically speaking, 
the lamellipodium has the topology of a ring lying between two curves, the 
outer one representing the leading edge and the inner one an artificially 
drawn boundary between the lamellipodium and the rest of the cell, roughly 
defined by a minimum actin density. More precisely, two non-identical inner 
boundaries for the two filament families are allowed. 

Actin filaments are polar with so called barbed and pointed ends. All 
barbed ends are assumed to meet the leading edge |[4Tj . 

Filament mechanics: Filaments are assumed to resist bending. More 
precisely, they are modeled as quasi-stationary Euler-Bernoulli beams. They 
are assumed to be inextensible [9]. 

Cross-links: The mechanical stability of the network largely relies on the 
existence of cross-links between the two families. Candidates for cross-linkers 
are proteins such as filamin [29], but also the Arp2/3 complex providing the 
connection between filaments at branch points [28] . It is assumed that cross- 
linking is dynamic with the building and breaking of cross-links as stochastic 
processes. While attached, cross-links are assumed as elastic, providing re¬ 
sistance against relative translational as well as rotational movement (away 
from an equilibrium angle) of the two filament families [43] , Characteristic 
life times of cross-links are assumed to be small compared to the dynamics 
induced by actin polymerization. The corresponding limiting process, which 
has been carried out in ED, leads to a friction model for the interaction 
between the filament families. 

Adhesion to the substrate: Transmembrane protein complexes with in- 
tegrins as their most important constituent provide adhesive connections be¬ 
tween the cytoskeleton and the substrate [ITS} I36j. Similarly to cross-links it 
is assumed that these adhesions are transient with relatively small recycling 
times, such that the averaged effect is friction between the actin network and 
the substrate. The short life time of adhesions is another questionable as¬ 
sumption, only satisfied for fast moving cells, where so called focal adhesions, 
i.e. large and very stable adhesion complexes, do not occur. 

The FBLM of [31] is still rather far from a complete description of all rel¬ 
evant processes. Some of these gaps are filled by the extensions below. Most 
importantly, the total number of filaments is kept fixed and their length dis¬ 
tribution is prescribed in m Here, filaments will be added by branching 
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and removed by capping and subsequent decomposition. The length distribu¬ 
tion will be determined by a quasi-equilibrium between polymerization and 
severing. In [31] cell size is regulated by a model for the effect of membrane 
tension. Here, this will be replaced by a contractive force in the cell center 
due to actin-myosin interaction. In certain applications it might be appropri¬ 
ate to combine these two effects. A further extension is a little speculative 
from a modeling point of view, but stabilizes the FBLM: We introduce a 
repulsive effect between parallel filaments of the same family, motivated by 
Coulomb interaction caused by the significant charges distributed along fil¬ 
aments. Finally, instead of a given fixed polymerization speed as in m, a 
model will be formulated involving both the effect of a chemotactic signal 
and of local leading edge bending. More details about these extensions are 
given in the following paragraphs. 

Polymerization and degradation: A desired polymerization speed is 
determined between a minimal and a maximal value, depending on the lo¬ 
cal concentration of an activator like PIP 3 , determined by a simple signal 
transduction model for a given chemoattractant distribution along the lead¬ 
ing edge. The desired polymerization speed is modified by the pushing force 
depending on the curvature of the leading edge (see Fig. ® 0 ' 

Several degradation processes of filaments are known. Aided by proteins 
of the ADF/cofilin family, they can depolymerize at the pointed ends El , or 
bigger pieces of actin can be removed by severing proteins like gelsolin [3j. 
We assume a severing process (see Fig. ER a mathematical description of 
which will lead to formulas for the filament length distribution, replacing the 
ad hoc approximations used in [31]. 

Branching and capping: In a lamellipodium, new filaments need to be 
created in order to maintain a polarized state. New filaments are nucle¬ 
ated by branching off existing filaments of the other family at or near to 
the leading edge (dentritic nucleation model [l28t [46] ). To form a branch the 
presence of the Actin-Related Protein-2/3 Complex (Arp2/3) at the mem¬ 
brane is necessary. Arp2/3 needs to be activated by nucleation promoting 
factors. Activated Arp2/3 is incorporated in the branches and later (e.g. 
upon filament degradation) reenters the cytoplasm, from which it is again 
recruited to the membrane [19] , Finally filaments can be capped at their 
branched ends by capping proteins [53], which blocks further polymerization 
(see Fig. [l]A.). The addition and removal of filaments had not been taken 
into account in [31]. 

Confinement: As a consequence of polymerization and adhesion, cells 
would spread indefinitely according to the model components described so 
far. In [31] cell confinement has been modeled as a consequence of membrane 
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tension. However, there is some experimental evidence m that confinement 
is mostly due to contractive effects caused by actin-myosin interaction in the 
rear of the lamellipodium. Myosin has the ability to utilize energy to move 
towards the barbed ends of actin and because of its bipolar structure this 
leads to contraction of antiparallel actin filaments. This mechanism plays 
an important role in cell movement, because it allows the cell to pull its rear 
parts [su Isa Em. 

We assume that the artificial inner boundary of the simulated lamel- 
lipodium is chosen such that the actin-myosin interaction takes place in the 
interior region not covered by the lamellipodium model. The lamellipodial 
actin network is assumed to be connected to an interior network of acto- 
myosin bundles providing a contractive effect (see Fig. eh 

Coulomb interaction: So far the model assumed no direct interaction be¬ 
tween filaments of the same family. Upon trying to understand the bundling 
of F-actin, it was discovered that F-actin, similarly to DNA, behaves like a 
polyelectrolyte [48l [49 , 50] . This means that F-actin is negatively charged 
(about 4 e/nm) at physiological conditions, hence there exists a repulsive 
force between the filaments. On the other hand, it has been shown that 
certain positively charged polycations, like divalent metal ions and basic 
polypeptides ]49], which act as counterions and neutralize the negative charges 
along the filament, promote filament bundling. 

As a modeling assumption, we introduce a repulsive effect between fil¬ 
aments of the same family (see Fig. [TJC) . The consequential inhibition of 
bundle formation in the lamellipodium is desired, since this is not our mod¬ 
eling goal at present (although it will be in future work). An additional 
motivation is the fact that a lack of coupling between filaments of the same 
family may lead to numerical instabilities, which can be avoided by the dif¬ 
fusive effect caused by repulsion. 
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Figure 1: New ingredients to the FBLM. 

The rest of this work is structured as follows. In the following Section 

1 the model of 1 T[] will be recalled. The new aspects will be introduced in 
Section [3j and the complete new model will be summarized in Section [4] 
Finally in Section [5] we will demonstrate the effect of the new terms and 
the potential of the full new model numerically. The power and flexibility 
of the model will become especially evident in simulations of the polariza¬ 
tion process induced by a chemotactic signal, of steady movement, and of 
a turning process. Movies of these simulations including visualizations of 
the stochastic filament dynamics are contained in the Supplementary Ma¬ 
terial. Examples of moving cell shapes, influenced by the response to the 
chemotactic signal, are shown. 

2 The Filament Based Lamellipodium Model 

Detailed derivations of the FBLM presented in this section can be found 
in |Til [51], first simulation results in [33] . and analytical results for the 
rotationally symmetric case in [32] , 

In the following, the superscripts + and — refer to the two filament 
families, also called clockwise and, respectively, anti-clockwise. These super¬ 
scripts will however be omitted, whenever we concentrate on one filament 
family. Quantities related to the other family will then be indicated by the 
superscript *. 

A semi-Lagrangian description of each family is used, where one coordi¬ 
nate a is a filament index varying on a torus (represented by a £ [0,2-/r)) 
because of the ring topology, and the other coordinate is the negative arc- 
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length s along a filament measured from the leading edge, i.e. s E [—L, 0] 
with the maximal (simulated) filament length L. Because of the inextensi¬ 
bility of filaments, s can also be seen as a material coordinate. However, this 
is only true for fixed time since, by polymerization at the leading edge with 
speed v(a,t) > 0 and the consequential inward flow of actin relative to the 
leading edge, a Lagrangian coordinate would be a = s + v(a,r) dr. For 
this reason the material derivative (i.e. the time derivative at fixed (a, a)) 

D t = d t - vd s 


will be used below. 

The main unknown in the model is F(a , s, t ) E M 2 which, for fixed a and 
t, represents a parametrization of the filament with index a at time t. By 
the inextensibility assumption 

| d s F(a, s, t)| = 1, for (a , s) E B := [0, 2ir) x [-L, 0], t > 0 , (1) 

it is an arc-length parametrization. More precisely, F represents the de¬ 
formation of all filaments in an infinitesimal a-interval, where the length 
distribution of the filaments is determined by the given function r/(a,s,t), 
whose value (between 0 and 1 ) is the fraction of filaments with length at 
least |s|. The assumption that all barbed ends meet the leading edge implies 
that r/(a, 0, t) = 1 and that r/ is monotonically increasing as a function of s. 
Note that rj(a, s, t) d(a, s ) can be interpreted as the total filament length in 
the infinitesimal coordinate volume d(a, s). 

At time t, the lamellipodium is represented by the set C(t) = £ + (i) U 
C~(t) with C±{t) = {F ± (o;, s,t) : (a, s) E B}. Note that C + (t ) and C~~{t) 
do note have to be identical. We request, however, that they share the leading 
edge, which can be motivated by the assumption of tethering of barbed ends 
to the membrane ("acto-clamping model", m 

{ F + (a,0,t ) : 0 < a < 27t| = {F _ (a,0,f) : 0 < a < 2i r} . (2) 

The artificial inner boundaries (s = —L) might be different. 

For the interaction between the two families, the points where filaments 
cross each other have to be described. This is done on the basis of two 
assumptions: First, there are no crossings of filaments of the same family, 
i.e. the map i ?=t (-,-,f) : B —> C±(t) is invertible. In particular, we assume 
that 


det(d a F ± ,d s F ± ) > 0 (3) 

holds, corresponding to a clockwise parametrization by a. Second, for each 
pair of filaments from different families there is at most one crossing, which 
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is transversal. We need representations of the set C + (t) n £ (t) in the 
coordinate domains. First we identify all pairs of crossing filaments: 

C(t) := {(a + , a~) G [0, 2ir) 2 : 3s ± (a + , a~,t) : 

F + (a + , s + (a + , a~, t), t) = F~(a ~, s _ (a + , a~, t), t)} . 

Then we define the parameter domains for both families such that 

F±(B±(t),t) = £+(t)n£-(t): 

B^{t) := { (a 1 * 1 , s ± (a + , a~, t)) : (a + , a~ ) € C(f)} C B . 

Note that, by our assumptions, the transformations 

(a + , a~) (->• (a 1 * 1 , s ± (a + , a~, t)) (4) 

from C(t) to Bq{£) are invertible and can be combined to transformations 
^ : B^(t) —>• B±(t) with the property 

F T = F^ o on B (for fixed t). (5) 

The positions and deformations of the filaments are computed on the basis of 
a quasi-stationary force balance obtained by minimizing a potential energy 
functional, which contains contributions from the bending of filaments, the 
stretching and twisting of cross-links, the stretching of substrate adhesions, 
and the membrane tension. This is coupled to age-structured population 
models for the distributions of cross-links and adhesions, assuming the build¬ 
ing and breaking of these connections as stochastic processes. The resulting 
model involves continuous delay terms, since, for the computation of the 
stretching forces, past deformations of the filaments are needed. With the 
model in this form, numerical simulations would be very expensive, partially 
also because it mixes different length scales. Whereas the effects of interest 
occur on the (/jrn) scale of the width of the lamellipdium, the stretching of 
cross-links and adhesions occurs on molecular (nm) scales. This implies that 
motion on the lamellipodium scale of the two filament families relative to 
each other and relative to the substrate is only possible, if the turnover of 
cross-links and adhesions is fast compared to other mechanisms (e.g. poly¬ 
merization and depolymerization). 

The corresponding limit has been carried out formally [31] and rigorously 
for a simplified model problem [22]. It leads to a friction model. The original 
idea seems to be more than 50 years old and has been formulated first for the 
derivation of models for rubber friction [42]. Recently it has been used for 
the modeling of the plastic reorganization of tissues due to cell-cell adhesion 
dynamics [38]. 

The assumption of fast turnover of substrate adhesions is reasonable for 
fast moving cells such as fish keratocytes, but certainly not satisfied for focal 
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adhesions, i.e. large stable and long lived adhesion complexes as found, e.g., 
in fibroblasts. The limiting procedure removes not only the delay terms, but 
also the coupling to the population models, which can be solved explicitly 
in the limit. 

The variational procedure involves the formulation of a Lagrangian func¬ 
tional where, besides the potential energy, also an account of the constraints 
© and ([2j is included. Its variation leads to a weak formulation of the 
problem for F: 


i-2ir 


0 = 


/r M (C - C 0 )_ 


d a F 

\d a F\ 


■ d a (6F) ^ Atether^ ' 


da 


J s=0 


+ / [,u B d 2 s F ■ dg (6F) + !i A D t F ■ 5F + A inext d s F • d s (6F)] r, d(a, s) 
J B 

+ [ [n s (D t F - D* t F*) • 8F 

Jc(t) 

=F//’(< / 9 - <p 0 )d s F ± ■ <9 S (<LF)1 ijq* d(a, a*), (6) 


for all variations 5F, where the first line contains contributions from the 
leading edge, and the convention = (— F 2 , F\) is used. The first 

term corresponds to the tension of the membrane with the total length of 
the leading edge, 


r2n /*27r 

C := \d a F + (a, 0, t)\ da = / \d a F~(a, 0, t)\ da , 

Jo Jo 

its prescribed equilibrium value Co, and an elasticity coefficient . The La¬ 
grange multiplier for the constraint (|2]) is a function defined along the leading 
edge denoted by Atethen an d v is the unit outward normal along the leading 
edge. The second line of ([6]) deals with forces within individual filaments: 
resistance against bending with bending modulus fi B , friction with the sub¬ 
strate caused by adhesions with a friction coefficient fj. A , and a tangential 
force due to the inextensibility constraint ([Tj) with Lagrange multiplier Ai nex t- 
The third and fourth lines of ([6j describe the effects of cross-links between 
the two families with a friction term caused by resistance against stretching 
of cross-links with friction coefficient fi s and a turning force term caused by 
resistance of twisting cross-links away from the equilibrium angle tpo. The 
angle between the filaments is determined by 


cos <£>(a + , a , t) = 

d s F + (a + , s + (a + , a - , t) ,t) ■ d s F~ (a~, s~ (a + , a - , f) , t) , 

and /j?" is the corresponding stiffness parameter. The above mentioned limit 
of fast adhesion and cross-link turnover provides explicit formulas for the 
coefficients / i A , /i s , and /r T in terms of mechanical and chemical properties 
of adhesion and cross-link molecules. 
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A strong formulation of the Euler-Lagrange equations requires to trans¬ 
form the domain of the last integral in ([6j to B. For this purpose, we use the 
maps from C{t ) to Bc(t) described above and introduce the modified friction 
coefficient and stiffness parameter 


/ m 

in B c {t) 

X o 

else, 


V T = 



da* | 
ds I 


in B c (t), 
else. 


(7) 


The strong formulation is then given by 

0 = p B d 2 s (rjd 2 F) + p A 7]D t F - d s ( rjX inext d s F ) 

+fi s r]rfAV ± d s (p T rj7)*{p - ip 0 )d s F ± ^ , 

where the computation of the relative velocity AV = Df F — D*F* o ip* and 
of the angle between the families, cos cp = d s F ■ ( d s F* o ip*), requires the 
transformation ip* between the coordinate domains. 

The corresponding boundary conditions are 


- p B d s ( r]d 2 F) + r]X inext d s F =F - <p 0 )d s F ± 

| 0 for s = — L, 

~ | TAtether^ + (C - C 0 ) + d a for s = 0, 

7] d 2 F = 0 for s = —L, 0. 


3 Modifications and Extensions 


Length Distribution and Filament Number Regulation 

In the model of im the filament number was conserved and the length 
distribution of filaments was prescribed with a fixed maximum length. In 
this section the model will be extended to include the effects of capping, 
branching, and severing on the filament number and length distribution. The 
results partially depend on the polymerization speed, the choice of which will 
be discussed below. 

The changes in filament numbers by branching and capping require a dif¬ 
ferent interpretation of the length distribution 7j(a, s, t ). For fixed s, 7](a, s, t ) 
will be considered as the number density of filaments of length at least — s 
in terms of a. Instead of the uniform distribution r/(a, 0, t) da = da of 
barbed ends, values of rj(a, 0, t) different from one are allowed. The density 
of barbed ends per leading edge length is then given by 


p(a,t) 


V 

\d a F(a, 0, t)\ 
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For the other family, the barbed end density p*(a* ,t) is defined analogously. 
With s = 0, the map between the coordinate domains (see ([5])) reduces to 
a map a*(a,t), and in the following p* means p*(a*(a,t),t). In the rest of 
this subsection we shall deal with fixed values of a. The dependence on a 
will therefore be suppressed for ease of reading. 

We start with the evolution of the number of barbed ends and assume 
that it depends on the barbed end densities per unit length: 

d t r/(0,t) = f(p,p*)\d a F(0,t)\ (8) 

where / (p, p*) is the change of barbed end number per unit length and time, 
modeling the effects of branching and capping at the barbed ends. Capped 
filaments become inactive and are assumed to be depolymerized very fast, 
such that they can be eliminated from the system immediately. 

It is instructive to rewrite ^ in terms of the length x = f \d a F\ da along 
the leading edge, instead of the Lagrangian variable a. With the lateral flow 
velocity Vi = f dt\d a F\ da (implicitly given as part of the filament dynamics), 
it can be written as 


d t p + d x {vip) = f. 


(9) 


Branching is assumed to be limited by the availability of activated Arp2/3 
complex at the leading edge with density a(t) (number/leading edge unit 
length). Its equilibrium value in the absence of branching is denoted by ao, 
the branching rate at equilibrium Arp2/3 density by Kb r , and the capping 
rate by n cap . This leads to the model 

, _ a -* 

J — ^br P ^cap P ■ 

ao 

With the rate c rec of recruitment and activation of Arp2/3 from the cyto¬ 
plasm, the evolution of Arp2/3 at the leading edge is governed by 


da 

dt 



-Kbr-(P + P*) , 
Oq 


The second term reflects the fact that activated Arp2/3 is incorporated into 
branches of both families. We assume that the Arp2/3 dynamics is fast 
compared to branching and capping and use the quasi steady (Michaelis- 
Menten) approximation 


®C)Crec 

Crec + Kbr (p + P*) 


f(p,P*) 


^br^rec P 

Crec + «br {p + P* ) 


^capP • 


The model (J9| has already been used in |TT] with a prescribed lateral flow 
velocity and with / = (5p*/(p + p*) — qp. A significant difference between 
our model and m is in the branching rate: For p = 0 and small values of 
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p*, the branching rate of HD is constant, i.e., not limited by the number 
of barbed ends of the other family, while in our model it is approximately 
n hr p*. 

An indication of the qualitative behavior of the model can be obtained 
from considering a situation, where the barbed end densities do not vary 
along the leading edge and are governed by the ODE system 

p = f(p,p*), p* = f{p*,p)- (io) 


It is easily seen that for Kbr < Av cap , he. capping exceeds branching, the 
densities converge to 0, otherwise the non-trivial steady state 


1 1 

K'cap ^br 

is stable. In separate work [21] we prove that this qualitative behavior carries 
over to the corresponding transport-reaction system with prescribed lateral 
flow velocities. 

The length distribution r/(s, t) of filaments is influenced by branching 
and capping through the model (|8]) for r/(0, t), but also directly by capping, 
which removes whole filaments by our above assumptions. We make the 
modeling decision that newly branched filaments are capped preferentially. 
This means that if branching exceeds capping, i.e., f(p,p*) > 0, we interpret 
/ as an effective net rate of production of new branches without further 
capping, whereas in the opposite case no new filaments are nucleated and 
—/is an effective rate of capping, affecting already existing filaments with 
a probability independent of their length. 

Degradation of filaments is also assumed to be facilitated by the action 
of severing proteins such as gelsolin [3] or ADF/cofilin m , cutting filaments 
at random positions. Similarly to the capping process, we assume that rear 
parts of filaments, which have been cut off, are completely decomposed im¬ 
mediately. 

Instead of using rj as unknown, it is more intuitive to write the model 
in terms of the density u(s, t ) = d s rj(s, t ) of filaments with respect to their 
length — s. Following mmm for the severing part, we use the model 




D t U = K sev / [H(s- s')v! - H(s' ~ s)«] ds' - Kcap ,eftU 


= 


«ds + SU - Atcap ,e&U , 


where v! abbreviates u(s',t). Cutting of a filament of length — s at the new 
length — s' occurs with the rate K, sev H(s r — s ) leading to a total severing 
rate —n sev s proportional to the length. This is the simplest possible model. 
More elaborate severing rates, e.g., with an account of aging of actin filament 
subunits ffl, could be included easily. 
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The effective capping rate is given, according to the discussion above, by 

(~f(p,P*))+ 

P 

It remains to rewrite the model in terms of r)(s, t) = u(a, t ) d a: 

D t rj = K S ev ST] K csl p t eSP ■ 

For a constant polymerization speed v this equation can be solved explicitly 
by the method of characteristics: 

r](s,t) = ??(0,t + s/n)exp -^ ^ K caPi e ff (r) dr^j , (12) 

where ?y(0, i) is given as the solution of ([ 8 ]). Assuming that the changes of 
the polymerization velocity are slow compared to the filament dynamics, this 
formula is a valid approximation also for time dependent velocities v(t). 

Finally, the maximum simulation length L of filaments is defined by a 
cut-off at small actin densities, i.e. a value ?/ m ; n for r]. With the rough ap¬ 
proximations of replacing rj(0,t + s/v) by r?( 0 ,f) and K C ap,eff(T) by K caPi eff(t), 
L can be computed explicitly: 

L(t) := - Kca P’ eff(t) + / Wff(^) 2 + lQg 7 ? (13) 

^sev y ^sev ^sev \ ^7min / 

resulting in the time dependent coordinate domain 

B(t) = {(a, s) : 0 < a < 2-7T, —L(a, t) < s < 0} . (14) 

Most notably, faster polymerization leads to wider lamellipodia. 

Myosin contraction 

We model the effect of interior contractile acto-myosin bundles as forces 
pulling the (artificial) pointed ends of lamellipodial actin filaments. Neglect¬ 
ing the effect of substrate adhesion of the interior region, the pulling forces 
are assumed to add up to zero by the action-reaction principle. 

The question of an appropriate direction for the pulling forces arises. We 
consider two scenarios: One where the contractile bundles pull the filaments 
tangentially, and one where they pull towards a central point, chosen as 
the center of actin mass (without having particularly strong arguments for 
this choice). Although one can argue for tangential pulling, which does not 
perturb the directional order of filaments in the lamellipodium, this choice 
has two disadvantages: If filaments get too tangential to the inner boundary 
of the lamellipodium, tangential pulling fails to control the size of the cell. 
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Secondly, tangential pulling is a slightly unstable process, since it might 
reinforce any small deflections of the pointed ends. For those reasons and 
in order to allow more flexibility of the model, we include a mixture of both 
choices. 

We again use the notation of dropping the superscript ± for the family 
under discussion and using * for the other one. The magnitude of the tan¬ 
gential force acting on the filament with index a is denoted by ftan(oi) and 
that of the centripetal force by /i n (a). We define 17(a) as the normalized 
vector pointing from the center of mass 

C M ■= J rjF d(s,a) (^j r) d(s,a) 

to the end F(a, — L(a )) of filament a. The forces can be included in the 
boundary conditions as 

~li B d s {jqd 2 s F) + rj\ inext d s F =F 7/??> r (<j6 - (j) 0 )d s F L 

= V (f tm d s F + ,f m V) , s = —L . 

We postulate a scalar positive quantity A, which measures the size of the 
contraction effect and which is chosen as A := fi IP (A c — Ao)+ with the area 
A c encircled by the lamellipodium and its equilibrium value Aq. 

The forces /tan(a) and /; n (a) are determined by the conditions that 

1. the total force should be close to the current contractility A. 



2 . it should be split between the tangential and centripetal contributions 
according to a weight 7 G [ 0 , 1 ], and 

3. the sum of all forces has to be zero. 

Mathematically this is realized by minimizing 


r 2n 


10 


T](S = —L) 


(/tan-7^ + (/in-(l-7M) 2 


7 1-7 

with the constraint 

r 2ir 

V (s = -L ) [f ta nd s F(s = -L ) + J\ n V] da = 0, 


da 


(15) 


/ 0 


giving 


/tan(a) = 7-4[ 1 -A t • d s F(a, -L(a))\ , / in (a) = {l-^)A[l~ti • 17(a)], 

(16) 


with 
I 1 = 


r2ir 


- 1 r2n 


7T \ r ZlT 

r] [7 d s F <g ) d s F + (1 — ^)V <S> V] da ] / 7 [7 d s F + (1 — 7 )!/] da . 

0 J Jo 

In Section [5] it will be shown that myosin pulling can effectively control cell 
size, and that the contraction force allows to produce moving cells. For these 
reasons we neglect the contribution of membrane tension. 
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Filament repulsion 

We consider a repulsive effect between parallel filaments caused by Coulomb 
interaction. The presence of mobile charge carriers in the cytoplasm leads to 
Debye screening with a typical Debye length in the range of a few nm, such 
that only local Coulomb interaction can be assumed, leading to a pressure¬ 
like repulsion effect. The electrostatic energy 


^pressure — J P*& d A 

is added to the potential energy, where C(t) is the area covered by the fila¬ 
ment family under discussion at time t (again dropping the superscript ±), 
p &A is the filament length in the infinitesimal area element dA, with 


V 

det (d a F,d s F) ’ 


(17) 


(where ([3| has been used) and $ is the electrostatic potential (see |39j for 
a derivation from a microscopic model based on individual filaments). A 
quasineutral approximation (justified by the relative smallness of the Debye 
length) and an equilibrium assumption for the mobile charge carriers result 
in a model for the electrostatic potential of the form <h = 3>(p). As example, 
the Boltzmann-Poisson model for the mobile carrier density leads to d> = 
^ P logGo), p P > 0. 

For the purpose of computing its variation, the electrostatic energy is 
written in terms of the quasi-Lagrangian coordinates, 


U, 


pressure 


[F\ = [ 

JB(t) 


<&(p)r) d(a,s), 


with the variation in the direction 5F, 


fiU pressure[F] 5F = [p(p) [det (d a F, 8 S (5F)) +det (d a (SF ), d s F)] d(a, s ). 

J B 

The considerations below will show that for stability reasons the pressure 
p(p) = &(p)p 2 has to be a nondecreasing function of the density p, which 
holds for the Boltzmann-Poisson model p(p) = p P p. Although we expect the 
pressure only to act in the direction orthogonal to the filaments, this con¬ 
sideration has not entered the discussion so far. However, the action of the 
pressure along filaments is eliminated by the (incompressibility) constraint 
10 - 

Model problem: It is instructive to look at a one-dimensional model prob¬ 
lem, where points with Lagrangian label a 6 M move along a line with 
positions x(a, t) E M (assumed a strictly increasing function of a). The den¬ 
sity of points is then given by p = l/d a x. The electrostatic energy takes 
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the form U pTessure [x] = f $(l/d a x)da. Its (L 2 -)gradient is given by d a p(p). 
If only the Coulomb interaction and friction (i.e. adhesions) are taken into 
account, the dynamics is governed by the gradient flow 

d t x = -d a p(p). 

With the continuity equation (in Eulerian coordinates) dtp + d x (pdtx ) = 0, 
this can be rewritten in Eulerian coordinates as 

dtp = d 2 x p{p) = d x (p(p)d x p). 

This is a nonlinear diffusion equation, where nonnegativity of the diffusivity 
p'(p) is necessary for stability. 

For the lamellipodium model, we may hope that the pressure term, by 
causing diffusion in the ct-direction, avoids intersections within a family and 
thereby stabilizes the system by ensuring that the modeling assumptions are 
not destroyed by the dynamics. This stabilizing effect is sometimes useful 
for numerical simulations as demonstrated in Section [5j 

Polymerization Rate 

Polymerization rates and, consequentially, polymerization speeds v(a,t) are 
subject to different regulatory mechanisms. We consider reaction to chemo- 
tactic signals, where the cell senses concentration gradients of a chemoat¬ 
tractant and translates them to varying polymerization rates. This leads 
to cell polarization and directed movement. The chemoattractant binds to 
receptors on the cell membrane that can trigger signaling pathways produc¬ 
ing intracellular gradients along the membrane reflecting the distribution of 
occupied receptors. For example, higher concentrations of PIP 3 have been 
observed towards chemotactic signals at the leading edge of Dictyostelmm 
discoideum and of neutrophils m- This in turn is expected to induce a local 
upregulation of actin polymerization mm- 

We consider constant planar chemoattractant gradients with the chemoat¬ 
tractant concentration S(x,y) = So + S±(x cos(</? ca ) + y sin(<£ ca )), where 5j 
determines the strength of the gradient and ip ca its direction. We model 
a chemotactic response independent of the strength of the chemoattractant 
gradient. A normalized internal quantity defined along the leading edge is 
given by 


S(F(a,0,t)) — min S(F(/3,0,t)) 

/3g[0,2tt) 

max S(F(/3, 0, t)) — min S(F(/3,0,t)) 

/3g[0,2t r) /3S[0,27r) 


Typically, PIP 3 is only observed at a part of the leading edge, possibly as 
consequence of a thresholding phenomenon of the signaling pathway. To 
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account for this, we choose a threshold c E [0,1] and define 


I(a, t ) 


dca(«,t)~C 

1 —C 

0 


for d ca {a,t) > c, 
else, 


( 18 ) 


which can be interpreted as a normalized PIP 3 concentration. The desired 
polymerization speed is chosen between prescribed minimal and maximal 
values: 


^opt(cb i) — ^min T 7(a, f) (u max T’min) • 

Finally, the polymerization speed is reduced by the force required to bend 
the leading edge outwards. On the other hand, due to filament tethering, we 
expect some acceleration of polymerization at leading edge segments which 
are curved inwards. These effects are described by an ad-hoc model for 
the polymerization speed v, depending on the signed local curvature n(a) 
(positive for convex leading edge regions): 

v = --- 

1 + exp {k/K re f) 


4 Full New Model 


The functions and variables used are the same as in the original model de¬ 
scribed in Section [2] To account for the different width of the lamellipodium 
for different filament indices a, we replace B by 

B(t) := {(a, s) |a E [0, 2ir), s E [—L(a, t), 0]} 


where L(a,t ) is given by ( |13| ), as already introduced in (14). 

The full weak formulation reads 

„ / 


0 = 


lB(t) 


H B d 2 s F ■ d 2 SF + fi A D t F • 6F + Ainext d s F • d s SF I r, d(a, s) 


IB{t) . 


y bending adhesion in-extensibility 

p(p) det (d a F, d s 5F) + det ( d a 5F , d s F ) d(a, s ) 


pressure 


/ 


+ 


IC(t) 


P S ( D t F - D* t F*)-5FTp T (</> - <k) d a F ± -d 3 SF pp* d(a, a*) 


y cross link stretching 

+ [ (ft,ndsF + f in V) • 5F V 


cross link twisting 


s=—L 


da 


Atether^ ' $F 




s=0 


da . 


myosin contraction 


tethering 
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The filament dependent magnitude of the inner pulling force /t an and /i n are 
given by (16) and p is defined by ©>• The filament densities r)(a, s, t ) are 
determined by (|8]) and (12). 


For the strong form, i.e. the Euler Lagrange Equations, the modified 
stiffness parameters are defined analogously as in ([7|. 

0 = p B d 2 s (? jd 2 F) - d s (ij\ine X td s F) + p A pD t F (19) 

+ d s (p^daF 1 -) - d a (p(p)d s F ± ^J 

± d s - Md s F^ + {D t F - D* t F '*) , 

where in the equation for F + , the derivatives of F~ are evaluated at 
(a - (a, s, t), s _ (o;, s, t )) and vice versa. The corresponding boundary condi¬ 
tions are 

-p B d s (pd 2 F) - p(p)d a F L + r?Ain ex t5 s F ^ rjr]*p T ((j) - <j> 0 )d s F ± (20) 

f V {ftan{a)d s F + fm{a)V(a)) for s = -L, 

1 TAtetherZ' for S = 0 , 

rjd 2 F = 0 for s = — L, 0 . 


5 Numerical Approach and Simulation Results 

In this section we sketch the numerical method for the solution of 
(20), described in more detail in 


(19), 


Simulations of model problems will 
demonstrate the effect of the new ingredients to the model introduced in this 
work. Finally, the full model is used to simulate the reaction of a cell to a 
chemotactic signal. 


Numerical Method 


The numerical approximation of (19), (20) is a formidable task. Filament 


families as described here are a new type of continuum, where both analytical 
and numerical approaches are still in early states of development (see [321 i33] 
for first results). The new numerical method used here will be described 
only briefly, since it is the subject of parallel work [20}, which we refer to for 
details. 

Instead of dealing with the time-dependent domains the equations 

for both families are transformed to the rectangular and time-independent 
computational domain [0, 27r] x [—1,0] by rescaling the variable s. A rectan¬ 
gular grid with uniform steplengths is used. The grid lines in the s-direction 
can be interpreted as computational filaments, each discretized by the same 
number of grid points, independent of its (time dependent) length. 
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The strong anisotropy in (19) is reflected in the choice of the finite element 
space for the spatial discretization. A tensor product space is used, where 
on each grid cell each component of F is represented by a fourth order 
polynomial, linear in a and cubic in s. These interpolate positions and first 
s-derivatives at the nodes. In other words, each computational filament is 
approximated by a cubic spline with linear interpolation in between. The 
finite element space is conforming for the weak formulation of (19), being 
continuous in a and C 1 in s and, thus, a subspace of H^{( 0 , 2i r), 1, 0 )). 

An implicit-explicit time discretization is used, based on a linearization. 
The evaluation of the interaction terms between the two filament families 
requires approximations of the mappings 'i/> ± , derived from Q, in order to 
represent filaments of one family on the grid of the other. The inextensibility 
constraint has been implemented by an Augmented Lagrangian approach. 


Effects of new model ingredients 

Filament Number Regulation: Branching and capping regulate the 
number of barbed ends. This is of particular importance when the poly¬ 
merization rate varies along the leading edge. For an explanation the lateral 
flow phenomenon is needed, i.e. the movement of barbed ends along the 
leading edge, caused by polymerization and the inclination of filaments. In 
our model, the (+)-family filaments typically move to the left (relative to the 
protrusion direction), and the (—)-family filaments to the right. Since bigger 
polymerization speeds increase the lateral flow, this would lead to filament 
depletion in regions with higher polymerization activity without the regu¬ 
latory effect of branching and capping. This typically reduces cell motility 
(see the paragraph on chemotaxis below and Figure [7]). 

To observe the regulation more directly, we started with a radially sym¬ 
metric cell, where the barbed end density along the leading edge is constant 
and equal to the equilibrium value ( |11[ ). Then the density of barbed ends of 
the (+)-family is perturbed locally. In the subsequent simulation a constant 
polymerization speed is used. Figure [2] shows the evolution of the barbed 
end densities for both families. If there is no branching and capping, the 
perturbation in r/ + is simply moved to the left by lateral flow without any 
changes (data not shown). The situation for ~p + is somewhat more compli¬ 
cated, because the geometry changes brought about by the higher number 
of filaments affects the number of filaments per length. Both ~p + and ~pT 
decrease initially because the membrane is locally pushed outward, making 
the cell slightly larger. The lateral flow is also visible here by the shift of the 
maximum filament number of p + . 

For the case where branching and capping are active, one can see how 
initially the number of (—)-family filaments increases because of branching. 
Additionally drops everywhere, again because the cell becomes slightly 
larger. However the dynamics eventually force the number of filaments to 
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return to its equilibrium value everywhere. 


Time:0.0 min Time:0.0 min 



Time:5.0 min Time:5.0 min 





Figure 2: Barbed end density perturbation with and without branching/capping: 
The density of left-moving filaments (red) is perturbed initially, that of right-moving 
filaments (dashed, blue) not. Thick lines represent the current state, thin lines 
represent the state at time t = 0. The left column shows the evolution in the 
absence of branching and capping, i.e. Kb r = K cap = 0. In the right column 
branching and capping are active with parameters as shown in Table [T] except 

H p = 0 . 


Actin-Myosin Interaction: Constraint (15) ensures that the myosin pulling 
on the inside of the lamellipodium is an internal force. We define Aq = r^ir. 
Figure [3] shows the evolution of the inner radius of a rotational symmet- 
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ric cell in time for different values of ?’o and p IP . It can be observed that 
smaller equilibrium areas and a stronger myosin force (i.e. larger p IP ) leads 
to smaller cell sizes. This shows that the actin-myosin interaction helps to 
control the cell size. In the section on chemotaxis below it will be shown 
that a moving cell can pull its rear due to this effect. 



Figure 3: Effect of the actin-myosin interaction: Comparison of the evolution of 
the radius of a circular cell for different values of ro and p IP . The horizontal lines 
mark the values of ro- All parameters as in Table [l] except the polymerization speed 
^opt 0.5. 


Pressure: The pressure term is a force acting only within one family. To 
demonstrate its effect we therefore look at a simplified model for one family 
only: We assume constant 77 , i.e. all filaments of the same length, no poly¬ 
merization or bending, and only the tangential component of the myosin 
force (7 = 1). Then the Euler-Lagrange equations simplify to 


H A d t F + d s (p(p)d a F ± ^j - 
-p(p)d a F ± + Ainext d s F 


da 


:F ± )- 


0 

ftan(a)d s F 


d s (A inext d s F) = 0 

for s = 0 
for s = —L(a) 


( 21 ) 


with the constraint |9 S P| = 1. To enable the computation of an analytical 
solution that can be compared to numerical results, we additionally assume 
the cell to be rotationally symmetric. This implies that the system can 
be fully described by a reference filament z(s,t). All other filaments are 
created by rotations. With the rotation matrix R(a ) for the rotation angle 
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a, F(a,s,t ) = R(a)z(s,t), leading to p = (z • d s z) 1 . Lastly we use the 
Boltzmann-Poisson model p(p) = p p p. For the pulling force we get / tan = 
p IP (irr 2 — 7T7 q) + . 

We introduce polar coordinates 2 = rexp(i</?) and denote by r and r' 
differentiation w.r.t. time and s respectively. This leads to the system 


p A r - p P “ A'nexl / - Ainext (r" - r(ip'f) = 0 . 

p A rip - A [ next ip'r - Ainext + <p"r) = 0 , 

(r ') 2 + (^')V = l, 

// P , _ f 0 for s = 0, 

' A mext 1 — S ,,IP 


Ainext ftp — / 


p IP (-Kr 2 — 7rrQ)_|_r' / for s = —L, 

0 for s = 0, 

p IP (nr 2 — nrfy+rip' for s = —L. 


( 22 ) 


The first two equations come from the Euler-Lagrange equations, the third 
is the in-extensibility constraint and the last lines stem from the boundary 
conditions. 

A stationary solution with straight radial filaments can be given explic¬ 
itly: 

p(s) = 0 , r(s ) = s + L + 77 , A(s) = p p ^1 + In 

(23) 


L + rj A 
s + L + ri) 


where the inner radius rj is the unique solution in (ro, oo) of 


ln (—“■) = » IP *( r l ~ r l) ■ 

As initial conditions in a simulation, we take a radially symmetric situation 
with straight, but not radial filaments. We then compare the numerical 
results to the analytical steady state given by (23). Figure |d] shows that 


convergence to the steady state is observed. This has been supported by a 
linearized stability analysis m Experiment with various initial conditions 
(data not shown) indicate that there might be an instability with respect 
to nonsymmetric perturbations. Since such an instability has never been 
observed in the full model, it was not investigated further. 
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Figure 4: Pressure. The upper pictures shows the evolution of the inner radius with 
time. The analytical steady state result ry (thick, dashed, red) and the parameter 
ro (thin, red) are also depicted. The stars mark times, where details are shown 
below: Here the filaments and Ai ne xt at three different times are shown. Again the 
analytical steady state solutions of the filament positions and of Ai nex t are shown 
(thick, dashed, red). [i B = fi T = /i s = 0, L = 10, A 0 = r = 157, /x p = 1. Other 
parameters as in Table [lj 


Chemotaxis 

Finally we want to demonstrate the potential of the full new model, in par¬ 
ticular its ability to simulate directed cell movement in the presence of a 
chemotactic signal, direction changes, and various cell shapes under differ¬ 
ent assumptions on the internal signalling network. We want to point out 
that we are aiming at a proof of principle and that serious comparisons with 
and fitting to experimental observations are the subject of ongoing work. In 
particular, although our parameter values (discussed below) and results are 
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in realistic ranges, we do not claim the simulated cell shapes to be close to 
experiments for any particular cell type. 


Starting to move: The following experiment mimics chemotaxis, i.e. a 
situation where the cell, in reaction to a chemical stimulus (chemoattractant) 
increases its polymerization rate (|12j.|54]). For the effect of the chemotactic 
signal on the polymerization speed we use the model introduced in Section 
[3] A time evolution starting with a rotationally symmetric cell put in a 
chemical gradient oriented to the right is shown in Figure [5] The first visible 
effect is that the lamellipodium on the right grows wider. The reason is 
that the increased polymerzation rate also increases the maximum filament 
length as modeled in (13). Next the cell starts moving, because filaments 
grow faster on the right. Additionally, the wider lamellipodium at the cell 
front exerts more friction than the thinner one at the cell rear. Eventually 
the cell shape remains constant and the cell moves steadily towards the right 
with a speed of about 3.9 fxm/min, which is in the biologically observed range 
(IS, 112]). In Figure [6] the two (numerical) steady states, the stationary cell 
and the moving cell, are shown together with some data. For the stationary 
cell the number of barbed ends is constant along the leading edge, whereas 
in the moving cell more filaments can be found at the back. To maintain 
this distribution the cell has to balance branching, capping and movement 
of filaments. For the stationary cell this simply means having branching and 
capping rates equal everywhere. For the moving cell, branching dominates at 
the front, whereas capping exceeds branching in the back. Additionally, the 
F-actin flow, i.e. the velocity of polymerized actin relative to the substrate 
is depicted in Figure [6] For the stationary cell, the flow is rather slow 
and uniform. For the moving cell one can observe small retrograde flow at 
the cell front and faster flow in the movement direction at the back, where 
retraction takes place. The flow speeds and distributions found are similar 
to those observed in HU 
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Time:0.00 min 



Time:5.50 min 



Time:11.50 min 



Time:15.50 min 



x(jifn) 
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Figure 5: Polarization and movement in the presence of a chemical gradient: A 
time series is shown, where the shading represents filament number, thin dashed 
lines show filament shape, and the thick filaments show the movement of a left 
moving (red) and right moving (blue) filament with time. Parameter values as 
in Table [T] with the internal signal threshold c = 0. For a movie of the simulation 
including a visualization of the stochastic filament dynamics see the Supplementary 
Material. 






















































































Barbed End Density F-Actin Flow 




5.5 

5 

4.5 


3.5 
3 

2.5 

2 

1.5 


5.5 
5 

4.5 

4 

3.5 

3 

2.5 
2 

1.5 


Figure 6: Nonmoving vs. moving steady state: Left side pictures show data along 
the leading edge. On the horizontal axes, 0 and 1 correspond to the cell front and 
0.5 to the cell rear. The upper left picture shows the barbed end density for the 
stationary steady state (blue) and the moving steady state (red). The lower left 
picture shows the branching (solid) and capping (dashed) rates. The two pictures 
on the right show the F-actin flow field, i.e. the velocity of polymerized actin 
relative to the substrate. Top: Stationary steady state. Bottom Moving steady 
state. Arrow length and color (values of the colorbars in fim/min ) represent speed. 
Parameters as in Table [lj 


Why filament number regulation? Figure [7] demonstrates the impor¬ 
tance of filament number regulation by branching and capping. The upper 
picture shows how for a cell with filament number regulation the filament 
densities at the rear and at the front remain close to each other and quite 
steady over time, whilst in the unregulated case they move away from each 
other. This is because in a cell where regulation is absent, filaments are 
transported to the back by lateral flow, which leads to an accumulation of 
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filaments at the cell rear, whereas the pulling front is depleted of filaments. 
The middle picture shows how this affects protrusion speed: In the regulated 
case rear retraction and front protrusion speeds approach the same value, as 
is necessary for constant movement, whilst in the unregulated case rear re¬ 
traction is slower and protrusion velocities decrease with time. In the time 
series below, one can see that this also affects cell shape: In the unregulated 
case filaments accumulate in the back, leading to a more prolonged shape. 
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Figure 7: Movement with and without filament number regulation: First row: 
Retraction speed at the rear (dashed) and protrusion speed (solid) at the font of 
the cell. Red lines (thick) show a cell, which can regulate filament number, blue 
(thin) lines a cell which cannot. Second row: Same as above, but showing the 
barbed end densities at the rear (dashed) and the front (solid) for the regulated 
(thick) and unregulated (thin) case. Third and fourth row of pictures: Cell shapes 
at different times. Shading represents actin density, thin gray lines the filament 
shape. Third row (blue, thin leading edge): Cell without regulation. Fourth row 
(red, thick leading edge): Cell with regulation. Parameters as in Table [l] except in 
the unregulated case Kb r = K cap = 0. 


A turning cell: One can ask the question if the steady state shape of a 
moving cell is affected by the initial conditions. A reasonable scenario for 
investigating this question, is a situation where the chemoattractant gradient 
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is gradually turned by 45 deg. Figure [8] shows the corresponding evolution. 
The final shape is very close to the original shape turned by 45 deg. 



Figure 8: A turning cell. The picture shows the cell shapes of a turning cell over 
9min. Top left inset: The direction of the chemoattractant gradient as a function 
of time. Bottom right inset: The initial state (red, dashed) and the turned final 
state (blue, solid) are compared. Parameter values as in Table [T] For a movie of 
the simulation including a visualization of the stochastic filament dynamics see the 
Supplementary Material. 


Various moving shapes: The shape of the moving cell strongly depends 
on the transduction of the chemotactic signal. In Figure [9]A-C three sce¬ 
narios are depicted, where only a certain fraction of the leading edge senses 
the stimulus. The more "local" the effect of the stimulus is, the longer the 
cell gets, because a smaller fraction of filaments pull the cell forward. The 
differences have been created in the model by variation of the threshold 
parameter c in (18). The fourth shape in Figure [9p shows results of an al¬ 


ternative mechanism, where not the polymerization speed but the branching 
rate is upregulated by the signal. The upregulating mechanism is as in [9]A. 
with a maximally threefold increase of the branching rate leading to a much 
denser actin network at the front. 


29 





























10 


Reaction to Chemoattractant 
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Time:31.66 min 









Figure 9: Different shapes of moving cells: Left pictures, A-C: Polymerization 
velocity along the leading edge (0 and 1 representing the cell front and 0.5 the cell 
rear). Left picture, D: Barbed end density along the leading edge. Right pictures: 
Final shapes. Shading represents actin density, thin lines in the cells (black, dotted) 
the filament shapes and lines at the leading edge (thick, blue) indicate the regions 
affected by the stimulus. A-C: The polymerization rate is affected, A: c = 0, B: 
c = 0.5, C: c = 0.7. D: The branching rate is affected, c = 0. Parameter values as 
in Table |T| 
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Parameters values: The parameter values used for the simulation are 
the ones summarized in Table [l] unless stated differently. Where availabe, 
values available in the literature have been used (e.g. for Aq, v±). 

Values for the branching and capping parameters Kbn K cap , c rec have been 
chosen in order to give 2p ref = 90 filaments per /im leading edge, which has 
been observed in real cells (|i!5]). K sev has been chosen in order to give a 
lamellipodial width in the relevant range of several fim [J5]. Other param¬ 
eters however (/i 5 , /i T , /i P , /i /p , K r ef) result from averaging processes in the 
derivation of the model (see [3Tj ). In principle they can be derived from 
molecular properties, but they depend in a complicated way on quantities, 
where not much experimental data is available, such as mechanical proper¬ 
ties of cross-linker molecules, their binding and unbinding rates, and their 
densities. These parameters can be used in a fitting process. Values within 
reasonable ranges have been chosen here. 
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Table 1: Parameter Values 


Var. 

Meaning 

Value 

Comment 


bending elasticity 

0.0 7pNpm 2 

m 


macroscopic friction 

caused by adhesions 

0.14 pNminpm~ 2 

measurements in |181 
[30|, estimation and cal¬ 
culations in [3H EH [3T] 

ftbr 

branching rate 

10 min~ l 

order of magnitude 
from |1TJ, chosen to fit 

2 p iei = 90 /im _1 [45] 

ftcap 

capping rate 

hmin~ l 

order of magnitude 
from [11 j , chosen to fit 
2/Vef = 90 pF 1 [45] 

Crec 

Arp2/3 recruitment 

900 pm~ 1 min~ 1 

chosen to fit 2 p ref = 
90 pm- 1 [45] 

ftsev 

severing rate 

0.38 min -1 /im,- 1 

chosen to give lamel- 
lipodium widths similar 
as described in |45| 

h IP 

actin-myosin interaction 
strength 

0.1pNpm~ 2 


Ao 

equilibrium inner area 

450 pm 2 

order of magnitude as in 

[521 M 

^min 

minimal polymerization 
speed 

1.5 pm/min -1 

in biological range 

%ax 

maximal polymeriza¬ 
tion speed 

8pm/ min -1 

in biological range 


pressure constant 

0.0 5pN pm 



cross-link stretching 

constant 

7.1 x 10 ~^pN min pm~ l 



cross-link twisting con¬ 
stant 

7.1 x 10 ~ 3 pm 


ftref 

reference leading edge 
curvature for polymer¬ 
ization speed reduction 

(5 pm)~ 1 
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